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Ontological theories of quantum mechanics describe a single system by means of well-defined 
classical variables and attribute the quantum uncertainties to our ignorance about the underlying 
reality represented by these variables. We consider the general class of ontological theories describing 
a quantum system by a set of variables with Markovian (either deterministic or stochastic) evolution. 
We provide the first proof that the number of continuous variables can not be smaller than 2N — 2, N 
being the Hilbert space dimension. Thus, any ontological Markovian theory of quantum mechanics 
requires a number of variables which grows exponentially with the physical size. This result is 
relevant also in the framework of quantum Monte Carlo methods. 



I. INTRODUCTION 

In classical mechanics, the number of variables describ- 
ing the state of a particle ensemble scales as the number 
of particles, thus the calculation time and memory re- 
sources which are necessary in a numerical simulation 
are generally a polynomial function of the physical size. 
This friendly property of classical systems seems to be 
lost at the microscopic level, when quantum effects be- 
come relevant. It is well-known that the definition of a 
quantum state requires a number of resources growing 
exponentially with the number of particles. This char- 
acteristic is at the basis of the exponential complexity 
of quantum mechanics (QM) and forbids at present to 
efficiently simulate the dynamics of many-body systems, 
unless some approximation is involved or particular prob- 
lems are considered [11, H, H, |1] ■ 

The origin of the exponential growth can be under- 
stood with the following example. A quantum scalar 
particle is associated with a complex wave- function ip{x) 
which lives in the three-dimensional space and evolves in 
accordance with the Schrodinger equation. Discretizing 
each coordinate with R points and evaluating the spatial 
derivatives by finite differences, the Schrodinger equation 
is then reduced to 2R^ real ordinary differential equa- 
tions. One could interpret ^ as a physical field which 
pilots the particle and intuitively conclude that M in- 
teracting particles are trivially described by M complex 
wave-functions in the three-dimensional space. Thus, the 
number of variables would grow linearly with the number 
of particles. Instead, in the standard description of quan- 
tum phenomena, an ensemble of M particles is associated 
with a single complex wave-function ip^xi, xm) which 
lives in the configuration space of the whole system. If 
each coordinate is discretized with R points, the number 
of real variables we have to integrate grows exponentially 
with M and is given by 2i?^*^. 

In general terms, if two quantum systems Ai and A2 
are associated with the Hilbert spaces Hi and 7^2, then 
the composite system Ai + Ai is associated with the ten- 
sorial product Jix^l-ii. This characteristic was at the ba- 
sis of the Born's argumentations against the realistic in- 
terpretation of the wave-function by Schrodinger, which 



would imply an exponential growth of physical variables. 
In the Born's interpretation the wave-function does not 
represent real quantities, but it is merely a mathemat- 
ical tool of the theory enabling to evaluate the proba- 
bilities of physical events. It is similar to the classical 
concept of probability distribution, that has the form 
p{xi,pi, ...,xm,Pm) for an ensemble of M particles and 
lives in the phase space of the whole system. In spite 
of this similarity, the Born's interpretation constitutes a 
deep departure from classical statistical theories, where 
each system is described by well-defined (ontological) 
quantities and the probability distributions merely pro- 
vide a statistical (epistemological) description of them (in 
the following, we use the terms "classical", "ontological" 
and "realistic" as synonyms). 

This departure is not unavoidable, indeed ontological 
theories of quantum phenomena exist, such as Bohm me- 
chanics d, la 0| I which describes single systems by means 
of well-defined variables (ontic variables) determining the 
outgoing values of measurements. However the Born's ar- 
gumentation constitutes the main difficulty of these ap- 
proaches, since they promote the abstract wave-function 
to the rank of a classical real field. Thus, it is nat- 
ural to ask whether there exists an ontological theory 
of quantum phenomena which uses an alternative rep- 
resentation of a single system and is not subject to the 
exponential growth of the number of variables. This is 
not only a foundational problem, but it has also a con- 
crete relevance, since such a theory would provide a new 
revolutionary quantum Monte Carlo method (QMC) for 
many-body problems [1, 0) 111- Questions pertaining to 
the computational complexity of QM C are discussed in 
Ref. f|. 

In this article, we consider this problem from a general 
point of view and give the first proof that the exponential 
growth of the number of variables is a common feature of 
any ontological Markovian (deterministic or probabilis- 
tic) theory of quantum systems. More precisely we prove 
that, for a system with a finite A^-dimensional Hilbert 
space, the number of continuous ontic variables can not 
be smaller than 2N — 2. As a consequence, the Born's 
criticism does not concern only the pilot-wave theories, 
but every realistic Markovian theory. If we assume QM 
exact, then this result suggests two choices, either to re- 
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ject ontological causal theories of quantum systems or to 
accept the exponential growth of the ontic variables as 
an intrinsic feature of quantum phenomena. After this 
work was completed, we became aware that similar ques- 
tions were discussed by L. Hardy 10]. He proved that 
the number of ontic states can not be finite, but did not 
placed constraints on the dimensionality of the ontic state 
space. As far as we know, the question on the exponen- 
tial growth of the ontological dimensionality was posed 
for the first time in Refs.jl, Q. Very recently, similar 
studies has been tackled in Ref. [llj. However, until now 
no answer to this open question was given. 

One could object that the use of the multi-particle 
wave-function as a variable field in ontological models is 
an obvious requirement, because of Bell's theorem on en- 
tangled states. Thus, our proof would not be necessary. 
However Bell's theorem imposes only a non-locality con- 
dition on realistic theories and says nothing about the 
dimensionality of the ontic state space. One can not re- 
ject a priori the possibility of a non-local realistic theory 
with a space of variables smaller than the Hilbert space 
without an explicit proof. 

In section |TT] we give the motivations for this study and 
introduce the general properties that an ontological the- 
ory of QM has to satisfy. In section iHll we provide some 
examples and in section IIVI discuss a particular case of 
dimensional reduction of the ontic state space. In sec- 
tion IV Al the proof of the theorem on the classical space 
dimension is presented and its consequences for quantum 
Monte Carlo methods are sketched. Finally, the conclu- 
sions are drawn in the last section. 



II. CLASSICAL THEORIES OF QUANTUM 
MECHANICS 



The state of a quantum system is described by a trace- 
one density operator p which is Hermitian and non- 
negative defined, i.e., it satisfies the properties 



Trp — 1 

{4'\p\4') ^ for every 



(1) 



When the maximal information on the system is ob- 
tained, the quantum state is described by a Hilbert 
space vector 1-0) and the density operator is the projec- 
tor A von Neumann measurement is associated 
with a Hermitian operator M = TOfe|fc)(fc|, where 
are the possible results and the vectors | k) are orthonor- 
mal. If M is not degenerate, when M is measured and 
the system is in the pure state |^), the probability of ob- 
taining result m/j is pk = \{k\ip)\'^. If is obtained, the 
quantum state is projected into |fc). This is the general 
framework of QM. 

In spite of the evident success of QM in explaining sec- 
ular experimental results, there are at least two reasons 
to ask for an alternative reformulation of the theory. The 



first one concerns the ambiguity in the state projection 
rule, which requires one to mark a boundary between 
the fuzzy microscopic quantum world and the macro- 
scopic well-defined observations. This ambiguity has at 
present no practical consequence, since the quantum pre- 
dictions are practically insensitive to the boundary, once 
the quantum domain is taken sufficiently large. This in- 
sensitivity is related to decoherence phenomena [12] . Our 
work is motivated by another more practical reason. In 
classical mechanics, the number of variables which spec- 
ify the state scales linearly with the physical size, thus a 
numerical simulation of dynamics is generally a polyno- 
mial complexity problem. Conversely, the definition of 
a quantum state requires an exponentially growing num- 
ber of resources, making the numerical integration of the 
Schrodinger equation impossible even for a small num- 
ber of particles. Many approximate methods are used 
in order to circumvent this problem, such as quantum 
Monte Carlo and semi-classical methods, the Hartrec- 
Fock approximation, the density-functional theory and so 
on. However, at present no general numerical method is 
known which is able to solve in polynomial time quantum 
many-body dynamics. It is interesting to observe that in 
the standard interpretation the wave-function does not 
represent a physical field, but it provides a complete sta- 
tistical information on an infinite ensemble of realiza- 
tions, thus the exponential growth of resources in the 
solution of the Schrodinger equation is not surprising. 
Indeed, this occurs also in classical mechanics with the 
definition of the multi-particle probability distribution 
p{xi^Pi, ...,xm,Pm)- For example, consider the prob- 
lem of M mutually coupled classical Brownian particles, 
which are described by the stochastic equations 



dvj 



+ ffi{t), 



(2) 



where Xi, Vi, Fij, 7 and ffi = {r]} ,r]'^ ,rif) are the spa- 
tial coordinates, the velocity, the interaction force, the 
coefficient of viscosity and the noise term, respectively. 
The two-time correlation function of the noise term is 
{ilf{t)rij{t')) = g5ijSki5{t — t'). The masses are set equal 
to 1. The solution of these equations is the trajectory 
of a single realization and its evaluation is a polynomial 
problem. Equation ([2]) is associated with the following 
Fokker-Planck equation 



dp 
dt 



E 



_d_ 

dvi 



d_ 

dxi 



(3) 

Its direct numerical integration is an exponential com- 
plexity problem. The difference of complexity between 
Eq. ([2]) and Eq. ([3]) is not amazing, since the first equa- 
tion describes the trajectory of a single reahzation, con- 
versely the second one provides a complete statistical 
description of an infinite number of realizations. Since 
this complete information is in general out of the experi- 
mental domain, a detailed evolution of the multi-particle 
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probability distribution is not required and it is practi- 
cally sufficient to evaluate the averages of some quan- 
tities over a finite number of trajectories by means of 
a Monte Carlo method. In QM , a similar approach is 
used for thermal equilibrium problems and, with some 
approximations, in dynamical problems (quantum Monte 
Carlo methods). The point of this discussion is the fol- 
lowing. Suppose, in accordance with the Einstein's view, 
that there exists a more fundamental theory which does 
not provide a statistical description on ensembles and 
characterizes each quantum system by means of a set 
of well-defined physical variables. The dynamical laws 
of this ontological theory and the Schrodinger equation 
would be similar to Eqs. (I2|3|) of our example, respec- 
tively. Thus, it is natural to pose the following question: 
is the trajectory simulation in this fundamental theory a 
polynomial problem? more precisely, does the ontologi- 
cal space dimension grow polynomially with the physical 
size? In this case, the fundamental theory would provide 
in a natural way a revolutionary quantum Monte Carlo 
method. This is the non-obvious core question of this 
article and we will find that the answer is negative, i.e., 
the exponential complexity is not related to the ensemble 
description, but it is a general feature of any ontological 
Markovian theory of quantum phenomena. In the fol- 
lowing subsections, we introduce the general properties 
of such theories. 



A. Kinematics 

We characterize a single system by means of a set 
of continuous and discrete ontological variables, say 
xi,X2,---^Xw and si,...,Sp, respectively. In the follow- 
ing we use synthetically the symbol X for this set, i.e., 
si, ...Sp}. For a single system, it takes a 
well-defined value X{t) at any time t. When a quantum 
system is prepared in a state the ontological variable 
takes a value X with a probability p{X) which depends 
on 1^). p{X) has to satisfy the following conditions 

J dXp{X) - 1 (4) 
p{X) > 0. (5) 

In principle, it is possible that does not fix unequivo- 
cally p; this one could then depend on the specific exper- 
imental setup used to prepare the pure state. Further- 
more in the quantum formalism the pure state prepara- 
tion deletes the memory of the previous history. This is 
not necessarily true in the ontological theory. In order 
not to lose in generality, we assume that each pure state 
can be associated with different probability distributions 
which can depend on the specific experimental setup and 
the previous history. This prevents us from writing a 
single- valued functional relation p{X\ip). Let C be 

a set which contains at least one element, we write the 



mapping 

|V)->MX|^,r7),77eC}. (6) 

This over-labeling of p can be found also in the posi- 
tive P- functions , mainly used in quantum optics and 
degenerate boson gases, and in Ref. where rj £ C 
is identified as context for the quantum state prepara- 
tion. Obviously, if \ip){i^\ ^ then p(X|i/',?7) ^ 
p(X\tp' ,7]') for any and 77' g C. Thus, Function © can 
be inverted and we have 

Wi^l^Dip), (7) 

where the operator D{p) is a function whose domain is 
the image of the functional p{X\4>,'i]). 

Equation ^ and its equivalent form, Eq. ([7]), are our 
first hypothesis. As a practical example, consider a sin- 
gle mode of the electromagnetic field, whose annihilation 
and creation operators are d and d^, respectively. The 
coherent state \a) of the mode is 

H=e-H^/^f:^|0), (8) 

where |0) is the vacuum state and a is a complex number. 
By means of the coherent state, it is possible to define 
some quasi-probability distributions associated with the 
quantum states. The Glauber distribution Paio:) of the 
state Itp) is defined as follows, 

j daPG{a)\a){a\ = \^){i^\. (9) 

A Glauber distribution exists for any quantum state, but 
in general it is highly singular and non-positive and can 
not be interpreted as a probability distribution [Eq. ([5|)]. 
The positive-P distribution P{a,P) is a generalization of 
Pg and is such that 

J dadf3P{a,(3)B{a,P) = \ip){^p\, (10) 

where B{a,P) = e(l"l'+l^l'-2/3*°)/2|a)(/3|. It is possible 
to prove that each quantum state \ip) is associated with a 
positive-P distribution [13,]. Equation (|10p is a concrete 
example of Eq. ([7]), X being the variable set {a, /?}. Note 
that a single mode has an infinite dimensional Hilbert 
space, conversely the {a, /?} space is four-dimensional, 
i.e, we can represent a quantum state as a statistical 
ensemble on a reduced space of variables. These vari- 
ables and P{a, (3) are analogous to the variables {xi, Vi} 
in Eq. ^ and the probability distribution p in Eq. ([2]), 
respectively. The dimensional reduction program would 
look, at this stage, feasible. However, we have still to 
define the dynamics of X and the connection between X 
and the measurements. 
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B. Dynamics 

In quantum mechanics, the state evolution of a conser- 
vative system from time Iq to ti is described by a unitary 
operator U{ti;tQ), i.e., we have the Markovian and de- 
terministic evolution 

mi)) ^u{h;to)mo)). (11) 

We retain in the ontological theory the Markovian prop- 
erty and define a conditional probability P{X,ti\X ,to) 
such that piiX) = J dXP{X,ti\X,to)po{X), where 
Po{X) and PiiX) are two probability distributions as- 
sociated with IV'(^o)) and \^{ti)), respectively. 

We assume that every unitary evolution is physically 
attainable. This hypothesis rests on the fact that, in 
quantum computers, every unitary evolution is in prin- 
ciple feasible by means of a finite number of quantum 
gates, that have a physical implementation (see for ex- 
ample Chapter 4 of Ref. [HI). As a consequence, every 
unitary operator has to be associated with a conditional 
probability. As for the probability distribution, in gen- 
eral each unitary operator can be mapped to many dif- 
ferent conditional probabilities which can depend on the 
the physical implementation of the evolution. Thus, we 
introduce a set £ which contains at least an element and 
define the mapping 

U-^{P{X\X,U,x),X(^£}- (12) 

The label yG £ identifies the context for the unitary 
evolution [14j . 

The conditional probabilities have to satisfy the fol- 
lowing conditions: 

1. For any U and x ^ ^ 

P{X\X,U,x)>0, (13) 
J dXP{X\X,U,x) = l. (14) 

2. For any [/, x G IV') and rj E C, there exists a 
rji E C such that 

p{X\Uij, m) = J dXP{X\X, U, x)p{X\i;, ry); (15) 

3. For any Ui^2 and xi.2 G £, there exists a X3 G ^ 
such that 

/ dX,P{X\X,, U2, X2)P{X,\X, C/i, xi) = .^g^ 

P{X\X,U2Ui,X3)- ^ ' 

Note that the integral symbol / dX synthetically indi- 
cates the integration and sum over the continuous and 
discrete variables. 

The set {P{X\X, C/,x),X G £} is an equivalence class 
labeled by U. The set of all the equivalence classes is a 
group, whose identity is {P{X\X, 1, x), X ^ 



These are the general hypotheses for the dynamics in 
the ontological theory. In the case of the positive-P dis- 
tribution, it is possible to define conditional probabilities 
which satisfy these conditions, but we will not show it. 

It is important for our purposes to deduce some prop- 
erties of P{X\X, U, x) and p{X\tl!, rj). We introduce the 
following 

Definition 1. We denote by I{ip,ri) the support of 
the probability distribution p(X\ip^ri), i.e., the smallest 
closed set with probability 1 and define lit/j) as follows 

/(7A) = {/(^,77),7yeC} (17) 

In practice, if the set of values X is countable, then X e 
I{'4','ri) ^ p{X\'ip,r]) ^ 0. For continuous spaces this is 
still true apart from a zero probability set. In order not 
to be pedantic, we assume that these negligible sets are 
null, i.e., we assume that our probability distributions 
have the same properties of discrete distributions. 

Property 1. li X E ^(V')j then the support of 
P{X\X,U,x) is a subset of I{U^). Equivalently, if 
X e I{ip) and PiX\X, U, x) 7^ 0, then X e I{Ui)). 

If the X space is discrete, this property is a conse- 
quence of Eqs (|5ll3ll5p . The proof is as follows: if 
X e /('f/') then there exists an -q such that p{X\il), rf) 7^ 0. 
Thus, if X is an element of the support of P{X\X, U, x), 
then exists an 771 (second enumerated conditions) such 
that piX\U^,m) = J:^P{X\X,U,x)p{X\^^,v) > 
P{X\X,U,x)piX\i^,r]) ^ 0, i.e., X is an element of 
I{U%jj). For continuous spaces. Property 1 is true apart 
from unimportant sets with zero probability. As previ- 
ously said, we will assume them null. 

Property 1 can formulated in these terms. If X{t) is 
the deterministic/stochastic trajectory of the ontologi- 
cal variable and X{tQ) G I{ip) at time then at a 
subsequent time ti > to X{ti) g I[U{ti;to)ip]- In gen- 
eral, this property is not invariant for time inversion and 
-'^(^o) i I{ip) does not implies that X{ti) ^ I[U{ti;to)ij]. 
The non-invariance is due to the fact that a backward 
Markovian process is not in general a Markovian process. 
Thus, the set I[ip{t)] is as a black hole, the trajectories 
X{t) can jump into it but can not escape from it. This 
is illustrated in Fig. [1^. The points (I), (II) and (III) in 
the figure are defined in the caption. 

It is reasonable to assume that this confiuence of the 
trajectories into I[U(ti]tQ)'ip] corresponds to a transient 
behavior which becomes negligible after a suitable series 
of unitary evolutions. This is equivalent to say that for 
any {ip'), there exists a probability distribution plXli/;' , rj) 
whose support contains only points of type (I) and (III) 
(Fig. [It). 

In the following, we assume that, after a suitable tran- 
sient evolution, the variable X evolves towards a clas- 
sical subspace where the time symmetry is fulfilled and 
consider only this subspace as the space of ontic states. 
Thus, we enunciate 

Property 2. For any \ip), U and x G if X ^ I{i/j) 
and P{X\X, U, x) 7^ 0, then X ^ I{Ui}}). 
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(a) . (b) 




FIG. 1: Visual representation of the evolution of I[^{t)] in 
the ontic state space, (a) Classical states in /[■i/)(fo)] at time 
to [labeled with (I)] can not escape from I[tp{t)] at subsequent 
time t > to. However, some states outside I[^p{to)] [labeled 
with (II)] can jump into I['4}{t)]. Other states (III) may exist 
which remain always outside I[il;{t)]. (b) The evolution of the 
support / for two different states and \^p). After a suitable 
transient evolution, no trajectory can escape from or to jump 
into I[i>{t)] and I[i>{t)]. 

Equivalently, if X{t) is a deterministic/stochastic tra- 
jectory and X(to) ^ liip) at time to, then at a subsequent 
time ii > to X{ti) ^ I[U{ti;to)^] (Fig. [lb). 

It is interesting to note that the actual state of a sys- 
tem we seek to describe in the laboratory might not have 
undergone the suitable transient evolution and not be in 
a region with the time symmetry property. For example, 
this might occur for a system sufficiently isolated since 
the early universe. However, it important to realize that 
the ontic state space is fixed once and for all and it must 
be able to describe every potential future evolution of 
the system. For our purpose, it is sufficient to know that 
there exists a region with the time symmetry property 
established by Property 2. The constraint on the dimen- 
sionality that we will find for this sub-region will be valid 
for the whole classical space. 

Properties 1-2 will be fundamental for our proof in 
Section |V3 

C. Measurements 

Let M be a Hermitian operator with eigenvectors \k) 
and eigenvalues mk, k being an integer index. When a 
measurement of M is performed on the state \4>), the 
result rrik is obtained with probability 

Pk = \{km' (18) 

if TOfe is not degenerate. 

In the ontological theory, we introduce a conditional 
probability Pjyj{k\X) such that 

5^P^^(fc|X) = l, (19) 

k 

PM{k\X)>0, (20) 



I dXP^,{k\X)piX\^p,rj) = |(fc|^)|2, (21) 

for any 77 G C and {ip). Equation ^T\\ makes the classical 
measurement rule equivalent to the Born's rule. 

When Pj^j{k\X) takes only the values and 1, the onto- 
logical theory is called " dispersion-free" and the variable 
X determines exactly the measurement results. However, 
this property is not necessary for our scope and will not 
be assumed. It is useful to note that each Hermitian 
operator M may be associated with many conditional 
probability distributions, i.e., the measurement results 
may depend on the physical implementation of the mea- 
surement of M. In order to account for this dependence, 
we have to introduce a set V, akin to C and and let 
the conditional probability depend on t € V, that is, de- 
note the conditional probability by Pj^-^(A;|X, r) [l3|. This 
function has to satisfy properties (|19l2ip for any t gV. 

For our purpose, it is sufficient to consider only trace- 
one projectors and prove the following 

Property 3. If j^') and \ip±) are two orthogonal 
states, then the supports of p(X|?/', 771) and p{X\'ip±,r]2) 
do not contain common elements, for every 771 and 772, 
i.e., n I{ip±) — 0. Proof: Assume ab absurdo the 
opposite and let Xo be a value such that p^Xolip, Vi) and 
p(A"o|7/'j_, 772) are not equal to zero. Let P(1|A") be the 
conditional probability of obtaining 1 with the projective 
measurement |7/')(7/'| and a fixed r. By Eq. (|2ip. we have 

J dXP{l\X)p{X\^,r,i) = l. (22) 

Since p{X\Tjj,rii) is positive and normalized to one 
[Eqs. (110], and p{Xo\ip,'ni) 7^ 0, by Eq. ^ we have 
that P(l|Aro) — 1. This is obvious if the set of values X 
is countable. As said in Section fll Bl this is true also in 
the continuous case apart from events with zero proba- 
bility which can be neglected. Similarly, we have 

J dXP{l\X)p{X\i,^, 772) = 0. (23) 

Since p(A"o|7/'j_, 772) 7^ 0, we have that P(l|Aro) = 0, in 
contradiction with the previous deduction. Thus, Xo can 
not be a common element of I{ip) and In simple 

words, since ji/)) and ji/jj,) correspond to mutually exclu- 
sive events, they can not be associated with the same 
value of the ontological variable. This property has been 
used in Ref. [l3| to derive a no-go theorem for noncontex- 
tual hidden variable models, and a similar property has 
been used in Ref. [10] to derive the "ontological excess 
baggage theorem" . 

It is interesting to note that every positive distribution 
introduced in quantum optics, such as the positive-P and 
the Husimi Q functions do not provide an ontological de- 
scription of quantum mechanics, since in general their 
support is the whole phase space and two orthogonal 
quantum states can have overlapping probability distri- 
butions. 
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We have completed our characterization of an ontolog- 
ical theory of QM. In the following section, we introduce 
some ontological models of simple quantum systems and 
show that their space dimension is always equal to or 
larger than 2N — 2, where N is the Hilbert space dimen- 
sion. 



III. TWO EXAMPLES OF ONTOLOGICAL 
THEORIES 

A. two-state quantum system 

We consider the ontological model of a two-state sys- 
tem reported in Ref. The classical variable X is a 
three-dimensional unit vector, which we denote by v. Let 
I — 1) and |1) be two orthogonal states, we associate each 
quantum state \ip) = "0-1 1 ^ 1) +"0111) with the following 
probability distribution in X, 

p{v\ij) = -V ■ w{tlj)e[v ■ w{-ip)], (24) 

TT 

where 9 is the Heaviside function and the three compo- 
nents of w{ip) are 

W2{i^) = -i^-iA + trii'-u (25) 
u;3(^)= |0-ip-|0i|2, 

i.e., w{ip) is the Bloch vector of w{ip) is the symme- 
try axis of the probability distribution, whose support is 
a hemisphere (the region where p is different from zero) . 
We write the Hamiltonian as 

3 

H^Y,hk{t)ak, (26) 

k=l 

CTk and hk (t) being the Pauli matrices and three real time- 
dependent coefficients, respectively. The vectors |±1) are 
the eigenstates of 0-3 with eigenvalues ±1. It is easy to 
prove that the equation of motion of w{ip) is 

^^2j2^^jkhjit)wk, (27) 

jk 

where e^fc is the Levi-Civita symbol. Thus, the probabil- 
ity distribution evolves rigidly according to the Louville 
equation 

^ = -^T.'^MtH^^, (28) 

ijk 

which corresponds to the deterministic equation for the 
ontological variable v 

^ = 2j2^rjkhjit)vk. (29) 



The measurement rule is as follows: the probability of an 
event associated with the vector \(p) is 

P{(j)\v) = 9[w{(t)) ■ v\. (30) 

It is possible to prove that 

J dvP{cf>\v)p{v\ij) = Mi;)\\ (31) 

the second side being the probability of the event 10) 
according to the Born rule. Thus, we can describe a 
two-state system as a classical one in a space with the 
dimension equal to 2N — 2, where = 2 is the Hilbert 
space dimension. 

It is useful to remark that, in this model, each vec- 
tor V is not associated with only one quantum state, 
i.e., probability distributions associated with different 
quantum states can overlap. In other words, if we de- 
fine S{v) as the set of Hilbert space vectors \tp) such 
that p{v\ip) ^ 0, then S{v) contains infinite elements. 
Note that this model satisfies property 3, i.e., if and 
\ip±) are two orthogonal vectors and jV') € ^{v), then 
\ip±) ^ S{v) (This property will be further discussed in 
Section |VT| . 

B. Higher dimension of the Hilbert space and 
Kochen-Specker theorem 

Classical dispersion-free models are possible also for 
higher dimensional quantum systems. Here we discuss a 
simple example introduced in the first chapter of Ref. [6| . 
It is very artificial, but shows that ontological formula- 
tions of quantum mechanics are possible. 

We consider a quantum system associated with an A^- 
dimensional Hilbert space. Let |1), |2),...,|iV) be an or- 
thonormal basis. We associate the quantum state 

(32) 

k 

with the probability distribution 

p{xi,--,XN,Mi^) = Y[5{Xk--ipk), (33) 

k 

where the ontic state space is spanned by the N complex 
variables Xk and A, which takes values in the real interval 
[0, 1] with uniform probability. Obviously, the ontologi- 
cal variables Xk evolves deterministically as ipk, i.e., by 
means of the Schrodinger equation 

^^^=E(^I^I0X^- (34) 

We can consider A a constant of motion. At this point, 
we have to write a conditional probability for events. 
Let \(j){l)),...\(j){N)) be a set of orthonormal vectors as- 
sociated with events. If the projective operator P^^^ = 
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(0(1)1 is measured, the probability of the event 



IS 



It is easy to prove that 



Pirn 

where 



d^^'x I dAP(l|x,A)p(x,A|V'), 





Pii\x..x)^o[\{x\<f>ii))\'-x] 



(35) 



(36) 



(37) 



with Ix) = J2kXk\k). Thus, the function P(l|x, A) can 
be interpreted as the conditional probability of the event 
1. If the projective operator \(j){2)) {(j){2)\ is subsequently 
measured, the probability of the event 0(2) is 



P(2|V>)^|(0(2)|V;)|^ 



(38) 



We want to find a conditional probability P(2|x, A) such 
that 

P{2\^) = J d'^x d\P{2\x, A)p(x, A|^). (39) 

Since the two events 0(1) and 0(2) are mutually exclu- 
sive, P(l|x, A) and P(2|x, A) cannot be different from 
zero for the same values of the conditional variables. 
Bearing this in mind, we put 

P(2|x, A) = 6[\{x\m)? + I(X|0(2))P - A] - P(l|x, A), 

(40) 

i.e., P(2|x,A) is 1 for |(x|0(l))P < A < |(x|0(l))P + 
|(x|0(2))P and zero elsewhere. Analogous construc- 
tions can be made for the other projective measurements 

\m){m\- 

It is interesting to note that the conditional probabil- 
ity of the event 0(2) depends by construction on |0(1)), 
i.e., a different choice of the first projective measurement 
modifies the outgoing result of the second one. This char- 
acteristic is called contextuality and is unavoidable when 
the Hilbert space dimension is higher than 2, as estab- 
lished by the Kochen-Specker theorem [iTj. 

We have shown that it is possible to construct an on- 
tological theory which fulfills the three conditions estab- 
lished in Section [III The classical variables are 2N + 1 
in number. However, since |xfeP = 1 f^nd the global 
phase is unimportant, the manifold of ontic states can 
be reduced to 2N — 1. Another dimension can be elimi- 
nated if we give up the dispersion-free property. A very 
simple example of model which is not dispersion-free and 
satisfies our three general conditions is obtained with the 
following probability distribution and conditional proba- 
bility for the state -0 and the event 0, 



p{xi 



(41) 



The corresponding ontological manifold has 2N — 2 di- 
mensions. Although this example sounds trivial, it shows 
that an ontological theory in a 2N — 2 dimensional space 
which satisfy our conditions is possible. As we will prove 
in Section IV Al this dimensional value is also the lowest 
possible one. 



IV. DIMENSIONAL REDUCTION OF THE 
ONTOLOGICAL SPACE IN A PARTICULAR 
CASE 

In this section we discuss an example of dimensional 
reduction of the ontological space. This reduction is pos- 
sible for particular quantum states and measurements. 
We will consider a bosonic mode and show that there 
exists a four-dimensional manifold in the Hilbert space 
whose elements can be represented in a two-dimensional 
classical space. Since this example is very well-known in 
literature, its discussion will be brief. For more details, 
see for example Ref. Q and references cited there in. 

We consider a one bosonic mode with a Hamiltonian 
quadratic in the annihilation and creation operators d 
and at. The results can be extended to the case of a 
higher number of modes. The Wigner distribution of a 
quantum state |^) is by definition the function 



W{a) = 



1 



d'\ (43) 



where a is a complex number and the domain of integra- 
tion is the complex plane. 

Let q and p be two Hermitian operators such that 
d = {q -\- ip)/\/2, they satisfy the canonical commuta- 
tion relation — i. In the basis of the eigenvectors 
I a;) of g. Equation (|43)) becomes 



W{q,p) r dxV(q + x/2)V'*(g-x/2)e-"P, (44) 



where q + ip = a and 0(x) = {xlijj). 

The Wigner function satisfies the identity 



J dqdpW{q,p) = 1. 



(45) 



In general, it can take negative values, but for particular 
states it is positive in the whole phase space (9,^) and 
can be interpreted as a probability distribution. This is 
the case of the Gaussian states 



■0(a;|go,Po,a, &) 



1 +ipo(x~qo)+ib(x-qof 



(7ra)i/4 



e 2" 



(46) 

whose Wigner distribution is the two-dimensional Gaus- 
sian function 

W{q,p\qo,Po.a,b) = ie-^^--[f-fo-2K9-?o)]^ (47) 



pm) = mx)?- 



(42) 



go and Pq are the mean values of q and p, a and b set the 
squeezing and the symmetry axes of the distribution. 
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From Eq. (|44|) it easy to verify that 

/oo 
dpW{q,p) = \^{q)\\ (48) 
-OO 

/OO 
dqW{q,p) = (49) 

where ■(/' is the Fourier transform of ijj. Thus, the 
marginal probabihty distributions of q and p are the 
probabiUty distributions of the observables x and p, re- 
spectively. In general, the probability distribution of 
the observables cos 9q + sin Op, called in quantum optics 
quadratures, is the marginal probability distribution of 
cos 9x + sin6'p. 

If only quadrature measurements of are considered, 
then the probability distribution of the outgoing values 
can be obtained by means of the classical probability 
rules of Section lll CI For example, we have from Eq. (|48|) 

\^{q)\^ = j dqdpP{q\q,p)W{q,p), (50) 

where P{q\q,p) = 6(q — q) is the conditional probability 
density associated with the measurement of q. Further- 
more, it is possible to prove that for any unitary evolution 
whose generator is quadratic in a and a^, the Wigncr 
function evolves according to the rules of Section IIIBI 
More precisely, the evolution equation of is a Liou- 
ville equation 

The manifold of the Gaussian quantum states is four- 
dimensional, conversely the number of classical variables 
is two. The dimensional reduction becomes more dras- 
tic when M modes are considered. In this case, the 
number of classical variables grows linearly with Af, 
conversely the dimension of the manifold of Gaussian 
quantum states grows quadratically. Multidimensional 
Wigner functions are used for example in the study of 
Bose-Einstein condensates [H, [l^ . 

V. THEOREM ON THE DIMENSION OF THE 
ONTOLOGICAL SPACE 

A. Definitions and properties 

In Section [n] we have established three general condi- 
tions of ontological theories of quantum mechanics, let us 
enumerate them. 

1. Let X, H and C be a set of discrete and/or contin- 
uous variables, a A''-dimensional Hilbert space and 
a set which contains at least one element. There 
exists a functional {p{X\tfj,ri),ri G C} which 
associates each quantum state lip) G H with a set 
of probability distributions of the variables X. 

2. Let £ he a, set with at least one element. There ex- 
ists a functional U — > {P{X\X, U, x), X ^ ^} which 
associates each unitary operator U oiTL with a set 



of conditional probabilities P{X\X ,U ,x)- These 
distributions satisfy the properties in Sect. Ill Bl 

3. Let P{(f)) be the trace-one projector \(t)){(f)\. There 
exists a conditional probability Pm{4>\X) such that 

j dxPM{(i>\x)p{x\i^,v) = {^\p{m) - i(0i^>i'- 

(51) 

for every rj and 

It is useful to introduce the following definitions: 
Definition 2. We call S{X) the set of quantum states 

IV') such that X e 
Definition 3. US{X) is the set of the vectors C/j-i/'), 

with IV') S S{X). In simple words, tJS{X) is the unitary 

evolution of the set S{X). 

By means of Properties 1-2, it is trivial to prove 
Property 4. If P{X\X,U,x) 7^ 0, then JJS{X) = 

S{X). 

Proof. It is sufficient to prove that a state ji/') is in 
S{X) if and only if U\^) G S{X). If |V') e S{X), then 
X e By Property 1, X e /(f/V'), i.e., J7|V') £ 

S{X). Similarly, it is proved by means of Property 2 
that U\il)) e S{X) ^ \^P) e S{X). 

Property 4 can be formulated in these terms. If X{t) 
is a deterministic/stochastic trajectory in the ontic state 
space and U(ti;to) is the associated unitary evolution 
from time to and to time ti, then 

Uituto)S[Xito)]^S[Xit,)], (52) 

for any io,i- An illustration of this property is reported 
in Fig. [21 where the unitary evolution U{ti]to) is syn- 
thetically denoted by U. At left we have represented a 
trajectory in a one-dimensional classical space. At right 
the corresponding evolution of S is sketched. 5[X(io,i)] 
are drawn as cones with I^E*) and [/|^') as symmetry axis. 
The Hilbert space is represented as a three-dimensional 
Euclidean space. 

Another fundamental property of S holds. 

Property 5. The set S{X) can not contain every 
vector of the Hilbert space. This is a consequence of 
Property 3. If \tp) is a vector of S{X), i.e. X G 
then any orthogonal vector of lip) is not an element of 
this set. 

At this point, the proof of our theorem on the classical 
space dimension is very simple. 

B. The theorem 

By means of the outlined properties, we will prove the 
following 

Theorem: The number of continuous variables in the 
set X can not be smaller than 2N — 2. 

In the proof, we will use Properties 4 and 5, which are 
a synthesis of Properties 1-3. 



9 




X space 



Hilbert space 



FIG. 2: Visual representation of a trajectory in a one- 
dimensional classical space and corresponding evolution of <S. 
At time to, the ontological variable takes the value X{to). At 
the subsequent time ti, the variable has evolved to a value 
X{ti) with probability P[X{ti)\X{to),U , x]- K this proba- 
bility is finite, the set S[X{ti)] of Hilbert space vectors (top 
right inset) is equal to the unitary evolution of S[X{to)] (bot- 
tom right inset), i.e., equal to US[X{tQ)]. The sets S are 
represented as cones with and as symmetry axis. 

Note that if the Hilbert space has a dimension higher than 2, 
the points of the one-dimensional classical space can not map 
every possible orientation of S[X{ti)]. 



Proof: Let X be a fixed value such that S{X) contains 
at least one vector of the Hilbert space. Property 4 says 
that for every unitary evolution U there exists a value X 
such that 



(53) 



This implies that the number of continuous ontological 
variables in X (xi, ...Xw; see Sec. Ill Al) is at least equal 
to the number of parameters required to specify the ori- 
entation of any unitary evolution of S{X) (see Fig. [5]). 
If S{X) would be invariant with respect to every U, one 
would find as lowest bound, but this is not our case, 
since S{X) is not a null set and does not contains ev- 
ery vector (Property 5). As discussed in the following, 
the lowest number of required orientation parameters is 
2N — 2. Thus, the theorem is proved. □ 

In an iV-dimensional Hilbert space the orientation of 
IJS{X) can be specified by iV — 1 orthogonal vectors, 
which evolve according to i.e., they are rigidly fixed 
in US{X). This representation is sufficient, but could 
be redundant when S has some symmetry. In order to 
clarify intuitively this point, one can consider some vi- 
sual examples in the Euclidean three-dimensional space. 
The orientation of a sphere does not require any param- 
eter, since it is invariant with respect to every rotation. 
A cylinder or a cone are invariant with respect to the 
rotation around their symmetry axis. Thus, in order to 
specify their orientation it is sufficient to give the direc- 
tion of this axis, i.e., a vector. By contrast, a pyramid 
has no rotational symmetry and we have to rigidly fix 



two orthogonal axes and specify their directions. In gen- 
eral, for a set of elements in an A'^-dimensional Euclidean 
or Hilbert space the orientation is specified by — 1 or- 
thonormal vectors when the set has no symmetry with 
respect to rotations, whereas no parameter is required if 
the set is completely symmetric. In our case, since S is 
not uniform, the number of required orientation vectors is 
at least one. This minimal labeling is possible when there 
exists one vector |\I') such that S is invariant with respect 
to every unitary evolution which leaves \'^) unchanged. 
In simple words, the orientation can be merely specified 
by one vector when there exists a symmetry "axis" and 
I^E") is the direction of this "axis" (see Fig. [2]). Thus, the 
set US{X) = Six) is identified by U\^), where l^-) is 
the symmetry axis of S{X). This vector is defined by 
N complex numbers, but because of the normalization 
of the quantum states and the irrelevance of their global 
phase, it can be labeled by 2N — 2 real parameters, which 
are the minimal number of parameters required to specify 
the orientation of US. 

Let us consider the two-state model of Section IlII Al as 
a practical illustration of the theorem. The space X is the 
set of three-dimensional unit vectors v and each quantum 
state is associated with the probability distribution in 
Eq. ([M)) . In this case, the elements oiS{v) are the Hilbert 
space vectors such that KV'lt/')^ > B = 1/2, |0) be- 
ing the state whose bloch vector, defined by Eq. (^5]) . 
is V. It is clear that \(f>) is the symmetry axis of S{v) 
and determines the orientation of the set. The dynam- 
ics of S{v) is a mere rotation and the axis evolves ac- 
cording with Eq. ([55)) . This model is dispersion- free. 
For more general models, it is possible to have differ- 
ent symmetric sets S{v) with < B < 1/2. However, 
B can not be larger than 1/2, because of Property 3. 
For B = 0, set S{v) contains only the vector |0) and the 
conditional probability defined by Eq. (|30|) is replaced by 
P(0|t7) = (1/2) [1 + iv{(j)) ■ v\. It is possible to have the- 
ories with an asymmetric set S, but in this case the two 
dimensions of the Bloch sphere are not sufficient. 

The hidden-variable model discussed by Aaronson in 
Ref. ^20] is another practical illustration of our result. 
Let |1),...,|A^) be a complete orthonormal basis of a iV- 
dimensional Hilbert space. In Ref. [20], the quantum 
state IV') = Q;i|1) -|- ... + unIN) is mapped to a probabil- 
ity distribution p{n) = a„ whose domain is a space of N 
ontic states. The unitary evolution is mapped to a, N xN 
stochastic matrix, which corresponds to the conditional 
probability P introduced in Section IIIBI As previous 
proven by Hardy [To], this mapping of quantum state 
does not work, since the number of ontic states has to be 
infinite. Furthermore, our theorem says that we need at 
least 2N — 2 continuous ontological variables. Aaronson 
makes up for such lack of a suitable number of ontic states 
by assuming that the stochastic matrix has to depend on 
the quantum state (See page 4 in Ref. [l^]). This cor- 
responds to assume that {ip) is an ontological variable, as 
in the Bohm theory. These models are typical examples 
of pilot-wave theories and the dimension of their ontic 
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state space is consistent with our constraint. 

Every known classical formulation of quantum me- 
chanics satisfies our constraint on the dimension of the 
classical space. In section IIIII and here, we have con- 
sidered some examples, other examples are provided in 
Refs. HH, S [13, [2I, where a iV-dimensional Hilbert 
space is reduced to a classical phase space of 2N real 
variables. Our result is relevant also for quantum Monte 
Carlo methods. By means of them, one tries to map 
a quantum dynamics to a classical one with a reduced 
phase space dimension. We have proved that this map- 
ping is not possible, unless some condition on the prob- 
ability distributions and the conditional probability dis- 
tributions is discarded. Quantum Monte Carlo methods 
in a reduced sampling space are introduced for example 
in Ref. [l|. In these methods, the kinematics and dy- 
namics conditions in Sec. |II]are fulfilled, but they do not 
provide positive conditional probability distributions for 
measurements. As a consequence, they are subject to 
the celebrated "sign problem" and have a computational 
complexity which grows exponentially with the evolution 
time and the number of particles. This complexity is not 
due to the dimension of the sampling space, but to the 
necessity of an exponentially large number of realizations 
in order to reduce the statistical errors. 



VI. CONCLUSIONS 

The proved theorem clearly shows that, if quantum 
mechanics is formulated as a Markovian realistic theory, 
the corresponding phase space grows exponentially with 
the physical size of the system. This occurs for theories 
which are able to describe every attainable unitary evo- 
lution. The explosion of the variable number implies that 
an exact Monte Carlo approach, which simulates quan- 
tum processes by means of realistic Markovian chains, is 
in general subject to an exponential growth of numerical 
resources and integration time [1^ . 

In some practical cases, polynomial algorithms could 
be feasible, for example when the trajectory of the quan- 
tum state is not dense in the Hilbert space or different 
states are indistinguishable experimentally. This is the 
case for example with dilute Bose-Einstein condensates 
in the mean field approximation. Decoherence could ac- 
tually play an important role in the complexity reduction 
for concrete systems. Roughly speaking, decoherence is 
due to our inability to distinguish a pure quantum state 
from a mixture of other quantum states ^26i] . When a 
system is composed of many particles, such as a macro- 
scopic gas, this inability is not merely technical, but fun- 
damental. Thus, a large set of states could be described 
as statistical mixtures of a smaller set of pure states, en- 
abling to simulate the dynamics by means of a Markovian 



theory in a reduced phase space, as expected in clas- 
sical limit. However, it is generally surmised that sys- 
tems as quantum computers require nearly unitary evo- 
lutions with negligible decoherence effects in order to ex- 
ecute efhciently quantum algorithms .ISi] and would not 
be suitably simulated by approximated stochastic meth- 
ods which replace pure state by mixtures. These con- 
clusions would support the conjecture that, in general, 
quantum algorithms cannot be efficiently simulated by 
classical computers in polynomial time [27] . Indeed, the 
actual speed-up of quantum computation requires fur- 
ther corroborations. For example, nobody has proven 
that factoring does not have a polynomial solution clas- 
sically. Note that the subset of available quantum states 
can be parametrized with a number of variables which 
grows with the number of quantum gates, i.e., with the 
physical size. If this subset is dense in the Hilbert space 
or its parametrization is computational hard, then our 
theorem supports the conjecture on the polynomial non- 
computability, but these conditions cannot be assumed a 
priori. The subset parametrization is certainly possible 
for example in particular classes of quantum computers, 
as recently reported in Refs. [2^[2^. 

Finally, we conclude with a discussion on some possi- 
ble implications of our theorem for the future study of 
the hidden variable theories. As said in the introduction, 
the main criticism against the known ontological mod- 
els, such as Bohm mechanics, is indeed the exponential 
growth of the ontic space dimension with the physical 
size. We have shown rigorously that this feature is un- 
avoidable in the framework of causal Markovian theories. 
Thus our result seems to put another nail in the hidden 
variable coffin and to imply that the realistic interpreta- 
tions do not provide a practical advantage for the study 
of quantum systems. However, our intent is to give a 
constructive result and to suggest a novel direction for 
the investigations on the ontological theories. In order to 
avoid the exponential growth of the number of ontic vari- 
ables, we have one possibility, to discard some hypotheses 
of the theorem. In our opinion, the Markovian property 
is the only one sacrifiable. More drastically, we could dis- 
card the causality hypothesis. It is interesting to observe 
that the Bell theorem and the Lorentz invariance seem 
to suggest the same conclusion. The Bell theorem estab- 
lishes that an ontological theory of quantum mechanics 
can not be local and relativity implies that a non-local 
theory is also non-causal. We suspect that realism in 
quantum mechanics with non-pathological consequences 
for the ontic space dimension could be possible with non- 
causal rules for evaluating correlations among events. 

I thank F. T. Arecchi for useful discussions. This work 
was supported by Ente Cassa di Risparmio di Firenze 
under the project "dinamiche cerebrali caotiche". 



[1] I. Camsotto, Y. Castin, J. Dalibard, Phys. Rev. A 63, 013618 (2006); R Deuar and R D. Drummond, Phys. 

023606 (2001); A. Montina, Y. Castin, Phys. Rev. A 73, 



11 



Rev. Lett. 98, 120402 (2007). 
[2] C. J. Trahan and R. E. Wyatt, J. Chem. Phys. 119, 7017 
(2203). 

[3] X. Oriols, Phys. Rev. Lett. 98, 066803 (2007). 

[4] A. Montina, J. Phys.: Conf. Ser. 67 012050 (2007). 

[5] D. Bohm, Phys. Rev. 85, 166 (1952). 

[6] J. S. Bell, Speakable and Unspeakable in Quantum Me- 
chanics, (Cambridge University Press, Cambridge, 1988). 

[7] D. Diirr, S. Goldstein, R. Tumulka, and N. Zanghi, Phys. 
Rev. Lett. 93, 090402 (2004). 

[8] A. Montina, Phys. Rev. Lett. 97, 180401 (2006). 

[9] M. Troyer and U-J Wiese, Phys. Rev. Lett. 94 170201 
(2005). 

[10] L. Hardy, Studies in History and Philosophy of Modern 
Physics 35, 267 (2004). 

[11] N. Harrigan, and T. Rudolph, larXiv:0709."TT49 [quant- 
ph]. 

[12] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 

[13] C. W. Gardiner and P. Zoller, Quantum Noise (Springer- 

Verlag, Berlin, 2000). 
[14] R. W. Spekkens, Phys. Rev. A 71, 052108 (2005). 
[15] A. Nielsen and I. L. Chuang, "Quantum Computation 

and Quantum Information" , Cambridge University Press 



(2001). 

[16] S. Kochen, E. P. Speaker, J. Math. Mech. 17, 59 (1967). 

[17] N. D. Mermin, Rev. Mod. Phys. 65, 803 (1993). 

[18] A. Montina and E. Arimondo, Eur. Phys. J. A 14, 105 

(2001); A. Montina, Phys. Rev. A 67, 043604 (2003); A. 

Montina, F. T. Arecchi, Phys. Rev. A 68, 053608 (2003); 
[19] M. J. Steel et al., Phys. Rev. A 58, 4824 (1998). 
[20] S. Aaronson, Phys. Rev. A 71, 032325 (2005). 
[21] F. Strocchi, Rev. Mod. Phys. 38, 36 (1966). 
[22] A. Heslot, Phys. Rev. D 31, 1341 (1985). 
[23] K. R. W. Jones, Phys. Rev. D 45, R2590 (1992). 
[24] D. J. Rowe, A. Ryman, G. Rosensteel, Phys. Rev. A 22, 

2362 (1980). 

[25] The exponential growth of the integration time is implied 
by the necessity of reading the whole information on X 
in order to evaluate the variables at the next time. 

[26] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 

[27] C. H. Bennett and D. P. DiVincenzo, Nature 404, 247 
(2000). 

[28] G. Vidal, Phys. Rev. Lett. 91, 147902 (2003). 
[29] N. Yoran and A. J. Short, Phys. Rev. Lett. 96, 170503 
(2006). 



